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MULTIPROJECTIVE WITNESS SETS AND A TRACE TEST 
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Abstract. In the field of numerical algebraic geometry, positive-dimensional solution sets of systems of 
polynomial equations are described by witness sets. In this paper, we define multiprojective witness sets 
which encode the multidegree information of an irreducible multiprojective variety. Our main results 
generalize the regeneration solving procedure, a trace test, and numerical irreducible decomposition to the 
multiprojective case. Examples are included to demonstrate this new approach. 
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1. Introduction 


Numerical algebraic geometry contains algorithms for computing and studying solution sets, called 
varieties, of systems of polynomial equations. Depending on the structure of the equations, we can view 
the variety as either affine or projective. Witness sets are the numerical algebraic geometric description 
of affine and projective varieties. If a variety is irreducible, its dimension and degree can be recovered 
directly from a witness set. When the variety is reducible, one can compute witness sets for each of the 
irreducible components thereby producing a numerical irreducible decomposition of the variety. 

We will consider multiprojective varieties, which are defined by a polynomial system consisting of 
multihomogeneous polynomials. Multiprojective varieties naturally arise in many applications includ¬ 
ing kinematics I41L likelihood geometry IITSl 1^ , and identifiability in tensor decomposition [14|. In 
fact, multihomogeneous homotopies f29l in numerical algebraic geometry developed from observing 
bihomogeneous structure of the inverse kinematics problem for 6R robots BTl . 

The paper is structured as follows. In Section [hi] we define multiprojective witness sets and their col¬ 
lections. Section summarizes homotopy continuation. Section contains our first main contribuhon: 
a membership test using multiprojechve witness sets. Section |4 contains our second main contribuhon: 
a generalization of the regeneration algorithm to compute multiprojective witness sets with examples 
presented in Section]^ Sectioncontains our third main contribution: a trace test for multiprojective 
varieties with examples presented in Section In particular, this trace test provided the motivation to 
recompute the BKK bound for Alt's problem (see Section [7^ in Theorem 7.1 


1.1. Multiprojective witness sets. A projective variety intersected with a general linear space has an 
expected number of solutions. When the dimension of the linear space is complementary to the variety, 
the expected number of solutions is finite and is called the degree of the variety. A witness point set for 
a variety is such a finite set of points. 

In the multiprojective setting, there are different types of linear spaces that may be taken which are 
related to the Chow ring (see (281 Chap. 8]). In particular, a collection of witness sets with all possible 
linear slicing can be used to describe a multiprojective variety. 
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Definition 1.1. Let V be a c-dimensional irreducible multiprojective variety in P”i x ■ ■ ■ x P”*" and let 
e = {ei,e 2 , ■ ■ ■ ,^k) S N>q such that c = |e| = ei + ■ ■ ■ + e^. An e* witness set for V is a triple: 

W® (V) := (V)} where 

(1) y is a set of polynomials that forms a witness system for V, i.e., V is an irreducible component of 

the solution set L = 0. 

k 

(2) L- = is a set of |e| linear polynomials with each being a general linear 

i=l 

polynomial in the unknowns associated with P”' of P”i x ■ ■ ■ x P”*^. The solution set L® = 0 
defines a codimension |e| linear space denoted by 

(3) w® (V) = V n £® is the witness point set for V with respect to L®. 

The number of points in w® (V), namely |w® (V) |, is the degree of V with respect to e, which we will 
formally denote by degW® (V) = |w®(V)|a;®. 


Definition 1.2. A (complete) witness set collection for a c-dimensional irreducible multiprojective variety 
V C P”i X • • • X P”*^ is a formal union of witness sets: 

2B(V)= □ W®(V) 

ecNio 

c=|e| 

with degree defined to be a formal sum: 

deg2B(V)= |w®(V)|a;®. 

e6N|o 

c=|e| 

When the context is clear, we will write W, W®, w® for 211 (V), W® (V), w® (V), respectively 


With this setup, the (multi)degree of V is deg2Il(V). 

Remark 1.3. One may disregard the terms where |w® (V) | = 0 in both the formal union of witness sets 
and the formal sum of degrees. 

Example 1.4. As an illustrative example, let V be the irreducible biprojective curve in P^ x P^ with 
coordinates ([xo,Xi], [i/O/i/i]) G P^ x P^ defined by V = {x\yo — x^yi). The degree of V and 211 is 
l(^(i,0) _|_ with the geometric meaning of this observed in Figure In particular, since V is an 

irreducible hypersurface, deg V can be observed directly from V based on the degree in each set of 
variables, i.e., deg V = deg,y{V)co^^''^^ + deg^{V)cv^^'^\ 


Notation 1.5. Our convention arises from a geometric interpretation based on slicing. This is the "recip¬ 
rocal" of the algebraic convention, which is followed by the multidegree function in Macaulay2 JTTJ. 
For brevity, the degree of a hypersurface Q computed by Macaulay2 is 


deg^ ^ |w 




,. . ., W 




Thus, the degree of the hypersurface V in Ex. 1.4 may be written as (2,1). 



MULTIPROJECTIVE WITNESS SETS AND A TRACE TEST 


3 


Figure 1. In the affine charf xq = 
I/O = 1/ the parabola intersects the 
vertical line at one point 

and the horizontal line at 

two points. 


yi 



Remark 1.6. According to Definition 1 1.2} w®(V) is a set of points. This set contains no information about 
multiplicity that is necessary for describing generically nonreduced components, i.e., components with 
multiplicity greater than one. One can attach the local multiplicity structure in the form of a Macaulay 
dual basis, e.g., ||9j[l2l[T3, to the points. Additionally, one can employ deflation methods, e.g., EOll^ . 
to perform computations on generically nonreduced components. 

Just as for classical witness sets, e.g., see jjSZl Chap. 13], we extend the above definitions to reducible 
varieties by taking formal unions over the irreducible components. 


2. Using Homotopy Continuation for Witness Sets 

A witness set provides information needed to perform geometric computations on varieties. The tool 
that permits such computations is homotopy continuation, which we briefly summarize in this section. 


2.1. Homotopies. Homotopy continuation is a fundamental tool in numerical algebraic geometry dis¬ 
cussed in detail in ||5l|37| and implemented in several software packages, e.g., lH|24l|^|38|. hr this 
manuscript, we employ straight-line homotopies, e.g., ||32l § 51]. One typical use of homotopy continu¬ 
ation below is to deform the linear space as in the witness set W®(V) to another linear space of the 
same type, say Ad®, along V. Let U be a witness system for V and suppose that £® and Ad® are defined 
by linear equations L® and M®, respectively. We denote this homotopy by 


H(V,£® ^ Ad®) 


V 

tU + {1- t)M^. 


The set of start points, at f = 1, for this homotopy is the witness point set w®(V). 

For membership testing (Section j^, one deforms to a general linear space of type e passing through 
a given point a*. The following specifies notation for such a linear space. 


Notation 2.1. Given a point a* G x 
passes through a*. 


X and e, let £®, be a general linear space of type e that 


Example 2.2. Let V C P^ x P^ and V as in Ex. 1.4 Let e = (0,1) and £® be the general linear space 


defined by L® = {yo — 5yi}. Then, w®(V) consists of two points, namely 

= (V) = {([v/5:l],[5:l]),([-v/5:l],[5:l])}. 


w 
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Consider cc* = ([1 : 1], [1:2]) and let which is defined by M® = L®, = {2yo — J/i}- The 

homotopy H(V, > TM®) = H(V,/I® —> >C®.) is thus 

1 =0 

t(yo - 5yi) + (1 - 0(2yo -yi) \ 


Starting, at t = 1, with the two points w®(V), the set of endpoints for this homotopy is 


{([l : V 2 ] ,[1 :2]),([l : -\/2] ,[1 :2 


For properly constructed homotopies, called complete homotopies in jlSl, each solution path defined 
by the homotopy, say Pi{t), is smooth on (0,1] and thus can be tracked using numerical path tracking 
methods, e.g., a predictor-corrector based approach. Endgames, e.g., see Il37l Ch. 10] and Il22l , are used 
to accurately compute the endpoints of the path, i.e., compute ^’^(O). 

Computationally, we perform path tracking in projective and multiprojective spaces by restricting to 
affine charts. That is, one imposes corresponding affine conditions on the coordinates thereby fixing a 
representation of the (multi)projective points. For example, in x P^, we can perform computations 
in C2 X C^ by taking affine charts of the form 


□xo -h Dxi -|- □ = 0 and Di/o + Clyi + 0 = 0. 

Here, we follow the convention of p5] where □ represents a random or unspecified complex number. 


2.2. Randomization. In a witness set for an irreducible variety V, the only condition on the polynomial 
system V is that it is a witness system (Def. |l.l[ ). As in Remark |l.6[ deflation techniques can be used to 
produce a system of pol}momial equations V such that V is a generically reduced irreducible component 
of V. That is, the dimension of the null space of the Jacobian matrix of V evaluated at a general point 
of V is equal to the dimension of V. In particular, the number of polynomials in V is greater than or 
equal to the codimension of V. When the number of polynomials in V is equal to the codimension of V, 
homotopies H(V, /I® —> Wl®) as constructed above are well-constrained systems, i.e., square. 

If the number of polynomials in V is strictly greater than the codimension of V, the homtopies are 
over-determined. As discussed in ||5j § 9.2], we will employ randomization for over-determined systems 
to improve numerical stability. That is, one performs numerical computations by replacing V with a 
generic randomization of V, say Randjy). Hence, V is a generically reduced irreducible component of 
Rand(y) in which the number of pol}momials in Rand(y) is equal to the codimension of V. 

In the multiprojective setting, we will maintain multihomogeneity by randomizing with respect to 
specified affine charts by fixing hyperplanes at infinity To that end, in each suppose that ELi is a 
general hyperplane defined by the linear polynomial H, . We construct the affine charts by having H, = 1 
and maintain multihomogeneity in the randomization via H,. This is demonstrated in the following. 


Example 2.3. Consider the following on x with variables [xo,Xi,X 2 ] and [yo/l/iW2]/ respectively: 

xfpi - xly2 

y ^ 

^?i/o - y^Q^iyi 
xIpq - XQXiyi 
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Let V be the solution set of V = 0, which is irreducible and has codimension 2. Hence, we aim to 
construct Rand(y) consisting of 2 polynomials by using the following: 

Hi = xq + 2xi — 3x2 and H 2 = 2j/o — 5yi + 3i/2- 

For example, we can take Rand(y) to have the form 

xfyi - xly2 + \JHi{xlyo - xoX2y2) + □Hi(x^yo - 


Rand(y)= 


xiX 2 yQ - XQyiy 2 + \JH 2 {xlyo - ^0X2^2) + C\H 2 {xlyo - XQXiyi) 
so that both polynomials in Rand(y) are homogeneous of degree (3,1) and (2,2), respectively 


One downside of utilizing randomization is the destruction of sparsity structure. For example, in 
Ex. 2.3 V consists of binomials while Rand(y) does not. Another downside of randomization is the 
increase of degrees. In or P^, one can order the polynomials based on degrees to minimize the 
degrees of the polynomials in the randomization, i.e., adding random linear combinations of smaller 
degree pol}momials to polynomials of higher degree. However, in the multiprojective case, the degrees 
of the polynomials, each of which is a vector of integers, need not have a well-ordering. 


3. Membership Test 


One application of witness sets is to use homotopy continuation to decide membership in the corre¬ 
sponding variety 13^ which was extended to images of algebraic sets using pseudowitness sets in [161. 
In this section, we describe our first main contribution which is using a multiprojective witness set 
collection to test membership in the corresponding multiprojective variety. 

Suppose that V C x • • • x P”* is an irreducible multiprojective variety and we want decide if a 
given point <x* is a member of V. The following highlights the difficulty of producing a membership 
test for multiprojective varieties, namely the loss of transversality of slices passing through a*. 

Example 3.1. Let V C P^ x P^ be the irreducible surface defined by 


y = 


xiyi - X2y2 
xoyiyi - xiyl 
XQy\ - X2yl 


and a* = ([1 : 0 : 0], [1 : 0 : 3]). Since V is the only irreducible component defined by V, we verify 
that a* G V by observing V{a*) = 0. 

The multiprojective witness set collection for V consists of three witness sets 

2B(V) = U U with deg2B(V) = + ico^o,2)_ 

In particular, for e = (2,0), the witness set W® has two witness points. For simplicity, let £® be 
the general codimension 2 linear space defined by L® = {—xq -|- 3xi — 2x2, Xq + xi + 3x2} and £®, be 
defined by L®, = {xi,X 2 }- The homotopy H(V,£® ^ >C®») defines two paths, both of which end at 
([1 : 0 : 0], [1 : 0 : 0]). Since a* is not an endpoint of this homotopy, a natural conclusion based on 
previous membership tests [Ml SSI is that a* ^ V. However, this perceived failure is obtained since 
these membership tests are based on the intersection of the variety and the linear space passing through 
the test point a* to be transverse at x* if a* is indeed contained in the variety. Here, V n >C®, is actually 
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a positive-dimensional set that contains a*, i.e., the intersection of the dimension 2 variety V and the 
codimension 2 linear space is not transverse at oc*. 

The following algorithm takes into account transversality for a membership test. 

Algorithm 3.2. [Membership test using multiprojective witness sets] Given a multiprojective witness 
set collection 2n(V) for an irreducible variety V C x • • • x determine if a given point <x* is 
contained in V. 

(1) For each W®(V) in 2If(V) with |w®(V)| > 0 

(a) Construct a general linear space of type e that passes through a*. 

(b) With start points w®(V), compute the endpoints E® defined by H(V,> /I®.). 

(c) If a* G E®, return "a* is a member of V." 

(d) If every point in E® is isolated in V H /I®, and a* ^ E®, return "a* is not a member of V." 

(2) Return "a* is not a member of V." 

Before proving correctness of Algorithm |3.2| we first show that, for each a* G V, there exists an 
element W®(V) G 2ir(V) with |w®(V) | > 0 such that a general slice >C®« is transverse to V at a*. 


Proposition 3.3 (Existence of transversal slices). IfV C P”i x ■ ■ ■ x P”* is an irreducible variety and cc* G V, 
there exists W®(V) G 2n(V) such that |w®(V)| > 0 in which a general linear space £®. of type e passing 
through CL* is transverse to V at cl*, i.e., cc* is an isolated point in V fl >C®.. 

Proof. For simplicity in our constructive proof, we write a* = {oc\,... ,0Lf) G P”' x ■ • • x P”*^ and set 

Mi — P”i X • • ■ X P”'-i X {a*} X P"-+i X ■ ■ ■ X 

Define Ci := dim V — dim(V fl A4i) and let /Ci C P”i x • • • x P”'' be a general linear space passing 

through CL* that imposes Ci conditions on P”i, i.e., /Ci is defined by Ci general linear equations in the 
unknowns of P”i that vanish on ccf We construct e 2 ,...,ei recursively. For i = 2,... ,k, define 

ei := dim (V n Wli n • • ■ n Wlf-i) - dim (V n Afi n • • • n Mi-i n Mi ), 

and let /Cj be a general linear space passing through cl* that imposes e, conditions on P”'. 

By construction, the intersection is a general linear space of type e that passes through cl*. 

Thus, we can take £®, to be It remains to show that for the constructed e, the linear space 

is transverse to V at a*. Since e, is the dimension of the fiber over cl* with respect to V fl /Ci fl • • • fl /C;_i, 
which is nonempty and has dimension at most n,, we know 0 < C; < n,. If c = dim V, then it 
follows that c = |e| because V is irreducible and a* G V. Since is sequentially imposing c; 

conditions on the corresponding fibers which have dimension e„ we know that cl* is an isolated point 

in V n /Cl n ■ ■ • n Upper semicontinuity, e.g., Ii37l Thm A.4.5], provides | w®(V) | > 1. □ 


The following illustrates this construction. 


Example 3.4. With the setup from Ex. |3.1[ we follow the proof of Prop. 3.3 First, 


Cl = dimV — dim (V Pi ({[1 : 0 : 0]} x P^)) =2 — 1 = 1. 

Hence, only one condition on the first P^ is needed to be imposed, say by /Ci as in Prop. 

62 = dim(V n /Cl) — dim (V n /Ci n (P^ x {[1 : 0 : 3]})) =1 — 0 = 1 


3.3 


Next, 
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showing that we also need to impose one condition on the second P^. Hence, e = (1,1) will yield slices 
that are transversal to V at a*. 

To illustrate, we consider the linear spaces and defined by 

L® = {xo+ xi + 3 x 2 ,i/o - -1/2} and L®. = {xi + 3 x 2 ,1/0 - - 1/2/3}, 

respectively The endpoints of the two solution paths defined by the homotopy H(V, —?■ >C®.) are a* 

and ([3 + 4-\/—1 : 3 : —1], [1 — 2\l —1 : 3]) with a* indeed being an isolated point of V fl 


We note that the constructive proof of Prop. 3.3 implicitly used an ordering of the spaces in the mul- 
tiprojective space P”i x • • • x P”''. If the spaces were ordered differently, one could yield a different e. 

Example 3.5. With the setup from Ex. 3.1 we follow the proof of Prop. 3.3 but take the P^'s is reverse 


order, i.e., constructing e = (61,62) by first computing 62 and then computing 61 . In this case, 

62 = dim y - dim (V n (P2 x {[1 : 0 : 3]})) =2-0 = 2 


and thus 61 = 0. Hence, e = (0,2) will also yield transversality at a* and therefore can also be used in 
a membership test. In fact, since w®(V) only consists of one point in this case, only one path is needed 
to be tracked to determine that a* G V. 


We use Prop. 3.3 to show that Algorithm |3.2|provides a membership test for multiprojective varieties. 


Theorem 3.6 (Correctness of Algorithm |3.2[ ). Algorithm |3.2| zs a valid membership test based on multiprojec¬ 
tive witness sets. 


Proof. If e is chosen such that a* G E®, then we know a* G V because E® C V H >C®, which justifies 
Item[H of Algorithm |3.2[ If every point in E® is isolated in V n >C®,, then coefficient-parameter homotopy 
theo ry |j30l provi des t hat E® = V n £®,. Hence, in this case, a* G V if and only if a* G E® which justifies 
Item[ld| By Prop. there exists e such that |w®(V)| > 1 and we have the property that a* G V if and 
only if a* G E® which justifies Item|^ 

□ 

Since Algorithm |3.2| requires the input of multihomogeneous witness set collection 2ir(V), the rest of 
the paper is devoted to computing such collections. 


4. Multiregeneration 

In this section, we will compute a witness set collection of a variety V. Thus, the aim is to compute 
isolated points of V n >C® for all possible e. When V is equidimensional, by all possible e we mean 
e G N>q such that |e| = dimV; otherwise, we mean 0 < e < {ni,...,nf). We take an equation- 
by-equation approach called regeneration IITSl [T9|. The key idea is the following: given a witness set 
collection for a variety y and a hypersurface Q, compute a witness set collection for y 'A Q. Iterating 
this process we will compute 2ir(V) with V defined as the intersection of hypersurfaces. Erom 21J(V), 
Section [ 6 ^ describes computing witness sets for the irreducible components of V. 

By utilizing an approach based on regeneration, one is actually constructing a sequence of witness 
point sets for subproblems as part of the computation. When the subproblems have physical meaning, 
the intermediate steps of regeneration provide useful information. Eor example, Alt's problem HI |Zl 
counts the number of four-bar linkages whose coupler curve passes through nine general points in 
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the plane (see Section 7.21. By using a regeneration-based approach, one actually solves the two-point, 
three-point, ..., and eight-point problems in the process of solving Alt's nine-point problem. 

The key step of multiregeneration is presented in Section 4.1 We summarize the complete algorithm 
in Section |4.2[ Examples are presented in Section 


4.1. Intersection with a hypersurface. Given a witness set collection for 3^ and a hypersurface Q, we 
compute a witness set collection for y r\Q. For the ease of exposition, we assume, in this subsection, 
that no irreducible component of y is contained in Q. The general case is provided in Section 

This computation has two steps: regenerate to union of hyperplanes and then deform to the hyper¬ 
surface Q. To simplify the presentation, we use the following. 


4.2 


Notation 4.1. We denote 5i to be the vector with 1 in the position and 0 elsewhere 


4.1.1. Regenerating to a union of hyperplanes. Following the notation above, let 3^ C P”' x ■ ■ ■ x and 

deg G = (^ 1 ,.. - rgk)- For i = 1,... ,k and ; = .,gi, let be a general linear pol}momial in the 

unknowns associated with P"'. Suppose that S is the union of hyerplanes defined by 

(1) s=nn»®- 

i=l ;=1 

Fet be the hyperplane defined by so that 

k Si ... 

s = u u 

With this setup, both G and S are hypersurfaces of degree {gi,. ■ ■ ,gk) and no irreducible component 
of y is contained in either G or S. From 231(3^), the following computes computes the isolated points 
in T” n 5 n for all possible d, thereby producing 221(3^ fl 5). 


Algorithm 4.2. [Regenerating to a union] Given a witness set collection 211(3^) and |sP| whose prod¬ 
uct S defines the union of hyperplanes S, compute a witness set collection 221(3^ fl 5). 

(1) Initialize 221(3^ fl S) to be the empty set. 

(2) For d E N>q such that there exist i E {1,2,... ,A:} and W^(y) E 221(3^) with d -|- = e: 

(a) Initialize = 0 and append to 221(3^ n 5) the set 

w‘^(3^n5) := |yu{s},l^p^}. 

(3) For i = 1,.. .,lc such that gi > 0: 

(a) For each W®(T’) G 221(3^) such that e > Sp 

(i) Define d := e — Si and hyperplane M. such that r\M. = CA. 

(ii) For ; = 1,.. . ,gi, append to the endpoints of the homotopy 


u{ync^,M^ 


'y 
< Ld 

tM+{l- t)s\’\ 
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with start points w®(3^), which are the isolated points of 3^ Pi 
(4) Return 211(3^ n 5). 


Remark 4.3. In Itemof Algorithm |4.2[ if 3^ is equidimensional then |d| = dim V — 1 and |e| = dim V. 
Proposition 4.4. Algorithm \4.2\ correctly returns the multiprojective witness set collection W{y fl 5). 

Proof. By construction, 

k gi ... 

3^n5n/:‘^ = ljlj3^n n 


In Item 


(/) 

3(a)i the hyperplanes A4 and Sf are both defined by general linear pol}momials in the un¬ 


knowns associated with P”c Hence, the isolated points in 3^ fl PI are obtained by the homotopy 


H (3^ P A4 —>■ sP ) starting at the isolated points of3^P£® = 3^PA4P 


□ 


4.1.2. Deforming hypersurfaces. The next step is to deform the hvpersurface 5 to ^ alone 3^ to compute 
2IJ(3 ^P^) from2B(3^P5). 

Algorithm 4.5. [Deforming hypersurface] Given 211(3^ P 5) and a hypersurface G such that S and G 
have the same degree and no irreducible component of 3^ is contained in G, compute a witness set 
collection 211(3^ P ^). 

(1) Initialize 211(3^ P to be the empty set. 

(2) For W^(3^ P 5) G 2B(3^ P S) do 

(a) Set Q'* to the set of endpoints of the homotopy H (3^ P C^,S G) starting at the isolated 
points of 3^ P 5 P 

(b) Remove from the nonisolated points of 3^ P ^ P e.g., via ||3l. 

(c) Append to 211(3^ P the multiprojective witness set 

W‘l(3^pg) := {yu{G},L^Q‘*}. 

(3) Return 211(3^ P^). 


Proposition 4.6. Algorithm \4.5\ correctly returns the multiprojective witness set collection W{y P G)- 

Proof. By Il30l , the set of endpoints of H (T* P C^, S ^ G) starting at the isolated points of 3^ P 5 P £** 
consists of a superset of the isolated points of 3^ P ^ P Hence, removing the nonisolated points, 
e.g., using |3|, yields the isolated points of 3^ P ^ P C^. □ 


Example 4.7. Consider V to be as in Ex. |l.4|with = {xq — 2xi} and = {yo ~ J/i}- In particular. 


2IJ(V) = {W(^'‘^)(V),W('^'^)(V)} where I (V) | = 1 and |w^°'^^(V)| = 2. Consider computing V P ^ 
where G is defined by 27xoi/o — SOy oXi — 25xoi/i -|- SOxiyi = 0 having degree {girgi) = (1/ !)• 

We can simplify Algorithm 4.2 by taking = xq — 2xi and = yo — y\. Hence, we have that 
2B(V P iS) = P 5)} where (V P 5) = (V) U w(‘^'^)(V) consisting of three points. 

For Algorithm |4.5} we only need to consider d = (0,0) for which Item 2a deforms from the three 
paths defined by deforming 5 to ^ along V as shown in Figure]^ The set of three endpoints is V P 
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Figure 2. Deforming from 5 fo 
along fhe parabola V in fhe affine 
charf defined by xq = i/o = 1- 


yi 



4.2. Multiregeneration algorithm. With Algorithms 4.2 and |4.5} we are ready to describe the full mul¬ 
tiregeneration process for a polynomial system V := |Gi,..., G^} which defines V C x • • • x 
The following computes witness set collections for in sequence. 

Algorithm 4.8. [Multiregeneration] 

Given V = {Gi,..., G^}, compute a wintess set collection 2I1(V). 

(1) Initialize a multiprojective witness set collection 2ir(P”i x ■ ■ • x P”*^) by solving a linear system 

of equations _ q 

(2) Fori = 0,...,£-1 

(a) Define A and Q := ^,+i. 

(b) Define y to be the union of irreducible components of X not contained in Q. 

(c) Define Z to be the union of irreducible components of X contained in Q. 

(d) Initialize the witness set collections 21J(3^) and W{Z) to each be the empty set. 

(e) For each W®(d:’) G W{X), do: 

(i) Partition w^{X) = N® U Lf® where N® = w®(d:’) \ Q and Lf® = w®(A’) fl Q. 

(ii) Append to 2B(i^) the witness set W®(3^) {X, L®, N®}. 

(iii) Append to 2B(Z) the witness set W®(Z) := {X U {G}, L®, Lf®}. 

(f) For deg^ = (^i,... ,gy), construct general linear forms Uf=i : / = 1,2,... and let S 
be their product as in 0 which defines the hypersurface S. 

(g) Input 2IJ(3^) and | into Algorithm 4.2 to compute 2IJ(3^ n 5). 

(h) Input 21J(3^ n iS) and Q into Algorithm 4.5 to compute 211(3^ fl Q). 

(i) Merge 2n(Z) and 2IJ(3^ n Q) to form 

(3) Return 2rr(V) = 

Proposition 4.9. Algorithm |4.8| correcf?y returns the multiprojective witness set collection 2n(V). 

Proof. For each e, we have N® = w®(3^) so that Algorithms |4.2| and |4.5| compute 21J(3^ fl ^). 

The set Z = X \ y is the union of irreducible components of X which are contained in Q. For each e, 
Lf® = w^{Z). Since each irreducible component of Z is an irreducible component of A n we know 
that each point in Lf® is an isolated point in A n ^ n G®. The result now follows immediately from 

xngnc^ = (yuz)ngnc^ = {yngnc^)u{znc^). 

□ 
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Remark 4.10. The performance of Algorithm |4.8| is based on the ordering of the polynomials and the 
structure of the polynomials. Example systems which can be viewed using different multihomogeniza¬ 
tions are considered in Section [S] 


Remark 4.11. As formulated, some of the paths which need to be tracked in Algorithm 4.8 may be 
singular. Following the notation in the proof of Prop. |4.9} this occurs when an irreducible component 
of y for i < £ has multiplicity greater than one with respect to {Gi,..., G,}. One can use deflation (see 
Remark |l. 6 | to produce paths with are nonsingular on (0,1]. 

Another option is to apply Algorithm |4.8| to a randomization of V (see Section In this case, all 
paths which need to be tracked are nonsingular on (0,1] by Bertini's Theorem so that one does not 


need to deflate paths. Two drawbacks of randomization, as discussed in Section 2.2 are the typical 


destruction of structure and the increase of degrees resulting in more paths which need to be tracked. 

The following is an illustrative example using Algorithm |4.8[ 

Example 4.12. Consider the following pol}momial system defined on x P^: 

V := {Gi, G 2 , Gs} := { xgyi - X2yi, xiy2 - X2yi, ^Qy\y2 - ^lyeyi } • 

It is easy to verify that V has four irreducible components: 

iSi = {{[xq,Xi,X 2 \, [1,0,0])} having dimension 2 and degree 
*52 = {([^O/^i/O]/ [j/O/J/i/O])} having dimension 2 and degree la;^^'^); 

= {{[xo,Xq,X 2 ], [i/ 0 / 1 / 0 / 1 / 2 ]) I X 2 yo = ^ 03 / 2 } having dimension 1 and degree Im^kO) + 

C 2 = {([ 0 / 0 / 1 ]/ [yo/ 0 , 1 / 2 ])} having dimension 1 and degree 
We demonstrate using Algorithm [4.8] to compute witness sets for the pure-dimensional components. 


i.e., 5i U S 2 and Ci U C 2 , with a decomposition computed in Ex. 6.5 

i = 0. Since A = P^ x P^ and ^ a hypersurface of degree (1,1), the witness point for P^ x P^ is 
regenerated into two points as summarized in the following chart. 


e 

(2/1) 

(1/2) 

|w®(^i)| 

1 

1 


L = 1. For both d = (2,1) and d = (1,2), = 0 so that all points are "nonsolutions" and must be 

regenerated. Since ^ = ^2 is a hypersurface of degree (1,1), two start points are regenerated into four 
points as summarized in the following chart. 


e 

(2/0) 

(1/1) 

(0/2) 

lw®(^i n ^ 2)1 

1 

2 

1 


L = 2. For d = (2,0), = 0 so that no points are regenerated. The point in is the witness point 

for iSi. For d = (1,1), both and consist of one point. The point in is a witness point for 1 S 2 . 
Since ^3 is a hypersurface of degree (1,2), regenerating this point yields two nonisolated endpoints, 
one each on iSi and 1 S 2 , and one isolated endpoint which is a witness point for Ci. For d = (0,2), 
consists of one point which must be regenerated. This yields two isolated endpoints, one each on Ci 
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and C 2 . Since V = H ^2 H ^ 3 , the results for V are summarized in the following chart. 


e 

(2,0) 

(1,1) 

(1,0) 

(0,1) 

|w'(V)| 

1 

1 

1 

2 


5. Multiregeneration Examples 

The following provides some examples demonstrating multiregeneration described in Algorithm |4.8| 

5.1. Comparison on a zero-dimensional variety. As mentioned in the introduction, multihomoge- 
neous homotopies 1 ^ were developed by observing a bihomogeneous structure when solving the 
inverse kinematics problem for 6 R robots. Following |j5l Ex. 5.2] with x denoting the cross product 
and o denoting the dot product, the polynomial system /(Z 2 ,Z 3 ,Z 4 ,Z 5 ) defined on (C^)^ is 

/i := zj o Z2 - Cl, /2 ;= Z5 o zg - C5 

ifS/ /s') (^ 1^1 ^ ^2 ^ 2^2 ^ 2^2 ^ ^3 ^ 3^3 ^ 3-^3 ^ ^4 “ 1 “ ^424 X Z5 -|- -|- U^Z^ X 2 g — p) 

ifeJjJs) ■= (Z2 0Z3-C2,Z3 0Z4-C3,Z4 0Z5-C4), if 9 , fw, fu, fu) ■= (Z2 o Z2 - 1 , Z4 o Z4 - 1 , Z3 o Z3 - 1 , Z5 o Z5 - 1) 

with variables Z 2 ,Z 3 ,Z 4 ,Z 5 G where zi,Z 6 , p G C^, ai,Ci,di G C are general constants. In particular, / 
consists of 12 polynomials in 12 variables and generically has 16 roots. We have selected the ordering 
so that the first two polynomials are linear, the next six are multilinear, and the last four are quadratic. 

We consider solving / = 0 using various multihomogeneous homotopies, a polyedral homotopy, 
and various multiregenerations. In particular, we consider using a 1-, 2-, and 4-homogeneous setup. 
Each setup corresponds with homogenizing / to yield a polynomial system V defined on a product 
of 1, 2, and 4 projective spaces based on natural groupings of the variables, namely {Z 2 ,Z 3 ,Z 4 ,Z 5 }, 
{z 2 ,Z 4 } * {z 3 ,Z 5 }, and {Z 2 } ★ {Z 3 } ★ {Z 4 } * {zs}, respectively. The corresponding multihomogeneous 
Bezout counts are 1024, 320, and 576, while the BKK root count is 288. These four counts are the total 
number of paths one needs to track using a 1-, 2-, and 4-homogeneous homotopy and a polyhedral 
homotopy, respectively. 

In the multiregeneration, we only use slices which could lead to isolated solutions yielding partial 
information about the multidegrees in the intermediate stages. Tables and summarize the multire¬ 
generation computations by listing the total number of start points and the total number of points in 
the resulting witness sets for each codimension. 

5.2. Comparison on a system with singular solutions. Our next example arises from computing the 
Lagrange points of a three-body system where the small body is assumed to exert a negligible gravi¬ 
tational force on the two other bodies, e.g.. Earth, Moon, and man-made satellite. The nondimension- 
alized system / defined on derived in ||51§5.5 . 1 ], with variables pi,w,di 3 ,d 23 ,x,y and parameter y, 
the ratio of the masses of the two large bodies, is 

f wpi - 1, WP2 - P, {pi - + y^- {P2 + xf -H - 5\y 1 

I (^<^13^23-‘^13)^ + 0F‘^23-P2^13' (®'^13‘^23 “ F'^23 “ ^ 13 )^ / 

where p 2 = 1 — pi- For generic y, this system has 64 solutions counting multiplicity: 32 of multiplicity 1, 
4 of multiplicity 2, 2 of multiplicity 3, and 2 of multiplicity 9. 

We first consider solving / = 0 using a multiregeneration applied to V, a homogenization of / based 
on using a 1-homogeneous and a 5-homogeneous structure defined by {pi} ★{^i 3 }*{ci 23 }* {x,y}. 
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1-homogeneous 1181 

2-homogeneous 

codim 

start points 

witness points 

start points 

witness points 

multidegree 

1 

1 

1 

1 

1 

la;{5.f>) 

2 

1 

1 

1 

1 

la;(5-5) 

3 

2 

2 

2 

2 

1^{4,5) + i^{5,4) 

4 

4 

4 

4 

4 

la;(3.5)+2a;<'‘''‘)+la;(5'3) 

5 

8 

8 

8 

8 

la;(2'5) +3a;(3'4) +3a;(4-3) + ^^( 5 . 2 ) 

6 

16 

16 

14 

14 

4u;P-4)+6u;(3'3)+4a;('‘'2) 

7 

32 

32 

20 

20 

10a;P,3) ^ i0a;(3'2) 

8 

64 

62 

20 

20 

20u;(3,2) 

9 

124 

90 

40 

34 

34 u;( 1 ' 2 ) 

10 

180 

62 

68 

28 

28a; (0'2) 

11 

124 

44 

56 

40 

40u;(0-i) 

12 

88 

16 

80 

16 

16a;(0'0) 

total 

644 


314 




Table 1. Summary of multiregeneration solving the inverse kinematics of 6R robot using 
a 1- and 2-homogeneous setup. 


codim 

4-homogeneous 

start points 

witness points 

multidegree 

1 

1 

1 

1^(2,3,3,3) 

2 

1 

1 

la;(2<3,3,2) 

3 

4 

4 

1^(1,3,3,2) i^(2,2,3,2) ^ i^(2,3,2,2) i^(2,3,3,l) 

4 

14 

11 

2^(1,2,3,2) ^ i^(l,3,2,2) ^ 1^(1,3,3,1) l^(2,l,3,2) 

+ 2a;(2A3'2) -f la;(3'2'3.i) + la;(2-3,i,2) + 2a;(2'3Ai) 

5 

28 

18 

3^(i,2,2,2) +2a;(iA3,i) + i^(i,3,1,2) +2a;(i'3-2,i) 

+ 3a;P'3A2) -f la;(3'i'3-i) + 3a;(2Ai'3) 

6 

18 

18 

6^(1,1,2,2) _^3^(1,1,3,1) +4^(1,2,1,2) +5^(1,2,2,1) 

7 

18 

16 

SwW.2) S^{1.1,2,1) 

8 

16 

14 

14^(1,1,1,1) 

9 

28 

24 

24a;(»'i'i-i) 

10 

48 

24 

24a;('>'0'i4) 

11 

48 

20 

20a;(0'0'0'i) 

12 

40 

16 

16^(o,o,o,o) 

total 

264 




Table 2. Summary of multiregeneration solving the inverse kinematics of 6R robot 
using a 4-homogeneous setup. 


This 5-homogeneous structure was selected since it minimizes the multrhomogeneous Bezout count 
amongst all possible partitions, namely 248 compared with the classical Bezout count of 1024. Even 
though the system has singular solutions, they only arise at the final stage of the multiregeneration 
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with this setup so that all paths which need to be tracked are nonsingular on (0,1]. The results are 
summarized in Table with the number of witness points counted with multiplicity. 



1-homogeneous 1181 

5-homogeneous 

codim 

start points 

witness points 

start points 

witness points 

multidegree 

1 

2 

2 

2 

2 

1^„{1,0,1,1,2) + 1^{0,1,1,1,2) 

2 

4 

1 

2 

1 

la,{0.04,l,2) 

3 

2 

2 

4 

4 

2^{0,0,1,1,1) + 2a;(0.0-0,i,2) 

4 

4 

4 

16 

14 

2^(o,o,i,i,o) ^ 4^(0,0,1,0,1) 4.4^(o,o,o,i,i) 4^{o,o,o,o,2) 

5 

32 

28 

48 

48 

40 a;(OAi,o,o) 40n;(0'0'»'i'«) + 28a;(0'»'0'0'i) 

6 

224 

64 

88 

64 

64^(o,o,o,o,o) 

total 

268 


160 




Table 3. Summary of using multiregeneration to solve the Lagrange points system 


Due to the presence of singular solutions, we also applied these two multiregenerations to the ho¬ 
mogenization of a perturbation of / as suggested in |2|, namely f + € where e G is general. From the 
isolated roots of / + e, the isolated roots of / are recovered using the parameter homotopy f + t ■ e = 0. 
The results matched those presented in Table with the only change being the endpoints in the final 
stage are all nonsingular. 


5.3. Comparison on a positive-dimensional variety. Our last example of multiregeneration arises from 
computing a rank-deficiency set of a skew-symmetric matrix. In particular, consider the system 










0 

Xl 

X2 

X3 

X4 

X5 









-xi 

0 

X6 

Xl 

Xg 

Xg 

S(x) 

•B- 

[ 1 

0 Ai 

A2 

A3 

A4 ■ 

n T 

where S(x) = 

-X 2 

-X 6 

0 

XV 3 

Xll 

X12 

S(x) 

■B- 

[ 0 

1 As 

Ae 

A7 

Ag ] 


-X3 

-X7 

-Xio 

0 

Xl 3 

Xl 4 









-X4 

-X 8 

-xn 

-Xl 3 

0 

Xl 5 









. -X5 

— Xg 

-Xu 

-Xl 4 

-Xl 5 

0 


and B G is a fixed general matrix. 

We focus on the skew-symmetric matrices S(v) of rank 4. That is, we aim to compute 


{(x*. A*) G X I rank S(x*) = 4 and A* is the unique solution to /(x*. A) = 0}. 

Since / consists of 12 polynomials and A G C®, the codimension of the projection onto x is at most 4. 
Therefore, from a witness set computational standpoint, we can intrinsically restrict to a 4-dimensional 
linear space in namely by fixing a general matrix A G and vector b G and considering 

the elements of of the form Ay -|- b where y G C"^. Thus, we consider the polynomial system 

F(y,A)=/(Ai/ + &,A). 

We consider polynomial systems Y arising from homogenizing F using a 1-homogeneous (C^^) and 
2-homogeneous (C^ x C®) structure. Since each polynomial in V, respectively, has degree 2 and multi¬ 
degree (1,1), we apply multiregeneration to a randomi z ation of V. In particular, the 1-homogeneous 
regeneration setup is equivalent to the regenerative cascade 1191 . Table ^summarizes the multiregen¬ 
erations for each codimension by listing the total number of start points, the total number of endpoints 
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codim 

1-homogeneous ll9l 

2-homogeneous 

start points 

iso. solns 

noniso. solns 

nonsolns 

start points 

iso. solns 

noniso. solns 

nonsolns 

1 

2 

0 

0 

2 

2 

0 

0 

2 

2 

4 

0 

0 

4 

4 

0 

0 

4 

3 

8 

0 

0 

8 

8 

0 

0 

8 

4 

16 

0 

0 

16 

16 

0 

0 

16 

5 

32 

0 

0 

31 

31 

0 

0 

31 

6 

62 

0 

0 

57 

57 

0 

0 

57 

7 

114 

0 

0 

99 

99 

0 

0 

99 

8 

198 

0 

0 

163 

163 

0 

0 

163 

9 

326 

45 

0 

178 

255 

45 

0 

178 

10 

356 

0 

84 

104 

260 

0 

84 

104 

11 

208 

0 

66 

24 

136 

0 

66 

24 

12 

48 

0 

20 

0 

24 

0 

20 

0 

total 

1374 




1055 





Table 4. Summary of various regeneration procedures for positive-dimensional solving 
related to the rank-deficiency system 


that satisfy the system and are isolated (U^), the total number of nonisolated solutions (removed in 
multiregeneration algorithm), and the total number of "nonsolutions" (N^). The only difference in 
this example between the 1- and 2-homogeneous regenerations are the number of start points of paths 
which need to be tracked yielding two additional costs for the 1-homogeneous setup. The first is the 
added cost in computing the additional start points themselves and the second is the typical added cost 
of tracking these extra paths which often become ill-conditioned near the end. 

In this example, the solution set of F = 0 is actually an irreducible component of codimension 9 and 
degree 45 in with multidegree -|- 24a;(^'^) in x C®. The closure of the 

projection of this irreducible component onto the y coordinates, namely 

{y G I rank S{Ay + b) < 4}, 

is a hypersurface of degree 3 that is defined by the cubic polynomial det S(Ay -|- b). 


6. Decomposition and a trace test 


After computing witness point sets for multihomogeneous varieties using multiregeneration de¬ 
scribed in Section]^ the last remaining piece is to decompose into multiprojective witness sets for each 
irreducible component. One key part of this decomposition is a trace test, first proposed in Il3^ for 
affine and projective varieties, which validates that a collection of witness points forms a witness point 
set for a union of irreducible components. We extend this to the multiprojective case in Section [O] The 
original motivation for such a test was to verify the generic number of solutions to a parameterized 
system with examples of this presented in Section 

Equipped with the membership test from Section th e trace test from Section |6.1[ and mon- 
odromy Il35l , we complete the decomposition in Section 6.2 


6.1. Trace test. Given a collection of witness points lying on a pure-dimensional multiprojective vari¬ 
ety V, the goal is to verify that they form a collection of witness point sets for some variety Z C V. 
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We accomplish this by generalizing the trace test proposed in Il36l for affine and projective varieties to 
the multiprojective setting. To that end, for V C x • • • x we first fix generic hyperplanes ; at 
infinity for each projective space P”' . Let H, be the polynomial defining T-Li, 

n n 

(2) H = YlHi and n = \Jni 

! = 1 i=l 

SO that "H is defined by the polynomial H. 

With this setup, we perform our trace using computations on the product space P”i x • • • x P”* with 
a general coordinate from the Segre product space p("i+i) '("*+i)-i defined as follows. 

Definition 6.1. A general coordinate in p("i+i) ' (")c+i)-i derived from P'^i x ■ ■ ■ x P”*^ has the form 

(3) P •= (n^o + • • • + • (Dyo + • • • + • • • (Dzo + • • • + Ozm^) 

where {x,y, . .. ,z) represents the coordinates in P”i x ■ ■ ■ x P”*. 

Hyperplanes in the Segre product space p("i+i) are defined by polynomials which are mul¬ 

tilinear, i.e., linear in the variables for each P”‘ . We say that TZ C P'^^ x • • • x P”*^ is a Segre linear slice 
defined by the polynomial R if R is multilinear. For example, in (|^ is a Segre linear slice since the 
polynomial H in (|^ is multilinear. With this, we are able to define the trace. 

Definition 6.2. Let V C P”i x • • • x P'^* be a pure c-dimensional multiprojective variety, 1 < c' < c, 
TZi,... ,TZci be general Segre linear slices defined by Ri ,... ,Rc', respectively, and be a linear space 
of codimension |e| = c — c' defined by L^. Let mCVn£®n7?.in---n TZc' and consider the homotopy 


(4) 


H® (VnL:^7^ ^ 7^ + 'H) 


'y 

< Ri + (1 - f) H 


l^R,,+ (l-f)H 

where H as in (j^. For each m G m, let m{t) denote the path defined by this homotopy starting at ni. 
Then, the trace of m with respect to H*' and a general coordinate p is the average of the p coordinate of 
the paths m{t), i.e., 

Trace^;,H®(Vni:^7^^7^ + 'H)(^) := ^ E P("^(0)- 

l“^l mem 

The trace is said to be affine linear with respect to "Hi,..., Tik if it is a linear function in t when restricted 
to the affine chart where Hi = • • ■ = Hjt = 1. 

We illustrate with the following example. 

Example 6.3. Consider the pure 1-dimensional variety V C P^ x P^ defined by 

^ = {(i/i - i/o)h - y\xo, {y\ - yl)x2 - y\xQ). 

In particular, V is the union of three irreducible varieties: 
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• Cl = {{[Q,Xi,X2], [1,1])}, 

• £2 = {([ 0 ,Xi,A 2 ], [1,-1])}, and 

• C = {{[xo,xi,xi],\yQ,yi\) \ {y\-yl)xi =y\xQ}. 

For simplicity, we take Hi = X 2 and H 2 = yi, and the Segre linear slice TZi be defined by 

Ri = (6Ao/7 + 3xi/5 + 2A2/7)(yo-i/i/2). 

The set V n 7?.i consists of five points nii,... ,m^, say where £1 n 7?.i = {mi}, £2 H 7?.i = {m 2 }, and 
Cr\TZi = {m 3 ,mi,ms}. With the general coordinate 

p - {2xo/7 - 5xi/12 + 3x2/17)(4yo/13 - 3yi/14), 


we consider the homotopy 


H® (V,7^l ^7^l + H) 


' {y\-yl)^2-ylxo 

_ ( 6 x 0 / 7 + 3xi/5 + 2 x 2 / 7 )(yo - 1 / 1 / 2 ) + (1 - t)x 2 yi. 


Restricting to X 2 = 1/1 = 1, Figure |^plots the trace of this homotopy for two sets of start points, namely 
{mi, m2, m3, mi} and {mi, m2, m3, mi, ms}. This plot shows that the trace of {mi, m2, m3, mi} is not affine 
linear while the trace {mi,m 2 , m 3 , mi, ms} is indeed affine linear. 



Figure 3. Plot of the trace for {mi, ..., mi} in blue (nonlinear) and {mi, ..., ms} in red (linear). 


As with the classical trace test l36l , the trace with start points m is affine linear if and only if m is a 
collection of witness points for a multiprojective variety. 

Theorem 6.4 (Trace test). Suppose that V C x • • • x is a pure c-dimensional variety defined by V, 
T-Li is a general hyperplane at infinity defined by Hi for P”‘, p is a general coordinate, TZi ,..., TZc' general 
Segre linear slices defined by Ri, .. .,Rc' for 1 < c' < c, is a linear space of codimension |e| = c — c' defined 
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hy L®, and m C V n mZiH ■ ■ ■ H TZc'. Then, Tracem H® (V n C^,7l TZ + TL) (t) is affine linear in t 
if and only if there exists Z C V n which is a union of irreducible components of dimension c', such that 
z n 7^l n ■ ■ ■ n 7^c/ = m. 

Proof. By applying Il36l Thm. 3.6] to the Segre product space ' one has the result if the 

trace in every coordinate of p("i+i) '(”t+i)-i is affine linear. Clearly, if the trace is affine linear in every 
coordinate of p("i+i) then it is affine linear in the general coordinate p. 

Conversely, we know, for all fia 2 r • • • / S C, 

(^ii ■ f^i2 ■ ■ ■ ~ ^ 

if and only if every = 0. Hence, if there exists 0, then 

■ ■ ■ Tk) Ihi-.k ^ ^ 

{k,-Ak) 

for general . ■ ■,Jak ^ C. Therefore, if the trace of some coordinate in the Segre product space 

p(wi+i)-•■(fjjj+i)-! jg gffij-^g linear in t, then the trace for p will also not be affine linear in t. □ 

In order to numerically test for affine linearity in t, a classical approach is to evaluate the trace 
at 3 distinct values of t and decide if they lie on a line. Another approach is described in @ which 
computes derivatives with respect to t of the trace. Two savings in this numerical computation are to 
work intrinsically on the hyperplanes H, = 1 to reduce the number of variables, and to use slices which 
preserve structure and/or simplify the computation. These savings are utilized in Section]^ 


6.2. Decomposition. Decomposition of a pure-dimensional variety corresponds with partitioning the 
collection of witness point sets into collections corresponding to each irreducible component. Similar 
to the classical trace test Il36l , one can use Theorem 16.41 to perform this partition by finding the smallest 


which has only 5 points to consider, this is an effective approach. However, this quickly becomes 
impractical as the number of points increases. Thus, we propose using two methods to reduce the 
number of possible partitions one needs to consider. Both are based on the fact that smooth points of 
irreducible components are path connected. 

The first approach is currently used in the classical affine or projective setting, namely to utilize 
random monodromy loops Il35l . That is, one tracks the points along the variety as a linear slice is 
moved in a general loop in the corresponding Grassmannian. Each start point and corresponding 
endpoint defined by this loop must lie on the same irreducible component. 

In the multiprojective setting, there is a second approach that one may utilize, namely to use one 
set of slices to perform membership testing (Section on the points arising from a different set of 
slices. As with monodromy loops, all points which are connected by smooth paths must lie on the 
same irreducible component thereby reducing the possible number of partitions to consider. 


subsets for which the trace is affine linear. When the total number of points is small, such as in Ex. 6.3 


Example 6.5. The multiregeneration computation in Ex. 4. 12| yielded pure-dimensional witness sets for 
dimensions 1 and 2. We now illustrate how to decompose into the four irreducible components. 
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We first start with the pure 2-dimensional variety which has degree Since the trace 

of m (2 0 )/ which is the unique point in V H is affine linear, this shows that there are two irreducible 
components, one of degree (namely, 5i) and one of degree (namely, 52 ). 

We now turn to the pure 1-dimensional variety which has degree For simplicity 

and concreteness, let Hi = X 2 and H 2 = yi define the hyperplanes as infinity Hi and H 2 , respectively 
Let = xq + x \ — 2 x 2 and = 1/0 “ ^l/i “ W define the linear slices and respectively, 
with Ri = define the Segre linear slice TZi. Thus, V n Hi consists of 3 isolated points 


nil = ([ 1 , 1 , 1 ], [ 1 , 1 , 1 ]), m 2 = i[-h-lA],[-h-lA]), and = ([ 0 , 0 , 1 ], [ 1 , 0 , 1 ]) 
where {wii} and {m2, m 3 } are the sets of isolated points of V n ^nd V fl respectively We now 
employ the membership test (see Section 31 applied to mi using For simplicity, we take = 


(1 -h V— 1 ) (1/0 ~ 2 i/i -|-1/2) which defines Cm[ and utilize the homotopy H(V, €,„[ ) with start 

points m2 and m3. The endpoints of the path starting at m2 and m3 are mi and ([ 0 , 0 , 1 ], [— 1 , 0 , 1 ]), 
respectively. Hence, mi and m 2 lie on the same irreducible component. The trace test shows that there 
are indeed two irreducible components, one from {mi, m 2 } having degree ImOT) -)_ (namely Ci) 

and the other from {m3} having degree (namely C2). 


7. Trace test examples 


We close with several examples using the multihomogeneous trace test from Section 


7.1. Tensor decomposition. In |8l, the problem of computing the number k of tensor decompositions 
of a general tensor of rank 8 in (8) C® ® is formulated with ||8[ Thm. 3.5] proving k > 6. Com¬ 
putation 4.2 of Ill4ll uses numerical algebraic geometry together with the fact the number of minimal 
decompositions of a general tensor of Sym^C^ (8 C^ (8 C^ is also k to conclude that k = 6. We use the 
trace test from Theorem |6.4| to confirm this result. To that end, we first consider a natural formulation 
as an affine variety. For a, & G C, Let ^ 3 ( 1 , a, h) G be the degree 3 Veronese embedding, i.e., 

Vo,{l,a,'b) = {l,a,b,a^,ab,b^,a^,a^b,ab^,b^). 

For X G C^, consider the tensor T(X) = ^ 3 ( 1 , Xi,X2) 8 ) (t3,X4) <8 ( 1 ,^ 5 ) G For general P, Q G 
we consider the set 


C = 


,xw. 


P.s + Q = J^T 


i=l 



C C X (C^) 


It follows from Il33l Thm. 3.42] that C is an irreducible variety of dimension 1. Using the natural 
embedding of C x (C^)^ x we have the two natural hyperplanes at infinity and can restrict 

each being set to 1 resulting in computations back on C C C x (C^)^. In order to maintain the invariance 
under the symmetric group on 8 elements, namely Sg, for ri,r 2 G C, we consider the Segre linear 
space Hi defined by 

Ri = (s - ri) H-^ 4 ^^ - ^' 2 ) . 

After randomly selecting P, Q G (C^)® and ri,r 2 G C, we used Bertini |4] to compute C fl Hi, which 
yielded a total of 528 • 8 ! = 21,288,960 points with precisely 6 ■ 8 ! points satisfying s — ri = 0, i.e., k = 6. 
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To verify that we have computed every point in C fl TZi, we applied the trace test from Theorem 6.4 
which showed that the trace was indeed affine linear confirming IIT4l Computation 4.2], 


7.2. Alt's problem. In 1923, Alt 111 formulated the problem of counting the number of coupler curves 
generated by four-bar linkages that pass through nine general points in the plane. The solution to this 
problem, namely 1442, was found using homotopy continuation in |40l. We use the trace test from 
Theorem 6.4 to confirm this result and study two related systems. 

Following the formulation in ||40| , only the displacements between one selected point and the other 
eight points are of interest. These displacements are represented using isotropic coordinates {Sj, dj) G C^ 
for y = 1,..., 8 . The polynomial system under consideration consists of 12 polynomials in 12 variables: 

a, b, m, n, x, y, a, b, m, n, x, y. 

The first four polynomials are: 

fi = n — ax, f 2 = n — ax, = m — by, = fn — by. 

Next, for; = 1,..., 8 , we have / 4 +j = 'yffj + + ^ 7 ° where: 

7 / 

Vj 



Tj 

- .yX .yj}/ y}/ yjX 

— ^jf^j 

^0 

■= and 


n - 5jX, 


n-5jX, 

J.X 

:= 5j[a — x) ^ 5j{a — x) 

3j5j, 

m — 3jy, 


:— m — Sjy, 

r] 

:= 5j{b-y)+Jj{b-y) 

- ¥;■• 


7 . 2 . 1 . Original problem. To study Alt's problem, we randomly selected {Sj,3j) G for j = 1,... ,7 and 
U(,V£ G C for ^ = 1,2. For = uis + Vi and = U 2 S + V 2 , we consider the irreducible variety 
C C C X defined by /i = ■ ■ ■ = /12 = 0 consisting of nondegenerate four-bar linkages whose 


coupler curve passes through {3j,Sj) for / = 1,...,8 . As in Section 7.1 we actually consider the natural 
embedding of C x ^ x and restrict each natural hyperplane at infinity to be 1. To further 
simplify the computation, we actually model the computation related to computing the closure of the 
image of C under the map 

n(s,a,b,m,n,x,y,a,b,m,n,x,y) = {s,x). 

That is, for random ri, r 2 G C, we consider intersecting C with the Segre linear space TZ-[ defined by 


Ri = (s-ri)(v-r2). 

We used Bert ini to compute C DlZi resulting in 32,358 points with precisely 8652 = 1442 ■ 2 ■ 3 of them 
satisfying s — ri = 0. In particular, 8652 corresponds to the number of distinct four-bar mechanisms 
while 1442 is the number of distinct coupler curves. The 6-fold reduction from mechanisms to coupler 
curves is due to a two-fold symmetry and Roberts cognates. Since the trace of the 32,358 points is 
indeed affine linear, this confirms the result of t40l. 
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7.2.2. Product decomposition bound. A product decomposition bound is constructed in 1311 by replacing 
for ; = 1 ,..., 8 with 

g4+j = {Oijlj + ) ■ imi + Vjl]) 

where ay, j 6 y, pj, vj G C. For generic choices, ll^ showed that the resulting system has 18,700 isolated so¬ 
lutions, which is a product decomposition bound on the number of isolated solutions for Alt's problem. 
We can repeat a similar computation for Alt's original problem with this product decomposition system 
by letting C denote the the union of the irreducible components of dimension 1 defined by randomly 
selecting all parameters except Ss- We then intersected C with the Segre linear slice TZi defined by 

Ri = {Ss - ri){x - r2) 

for random ri,r 2 G C. Using Bertini, we find that C HlZi consists of 37,177 points, precisely 18,700 
satisfy Ss — ri= 0 . We confirm the result of USTII since the trace of the 37,177 points is affine linear. 

7.2.3. Polyhedral bound. In our last example, we consider the family of systems iF with the same mono¬ 
mial support as the polynomial system fi ,.. .,/i 2 defining Alt's problem as above. The number of 
solutions to a generic member of iF is called the polyhedral bound or BKK bound Ii39l . We consider a 
line in T by randomly fixing all coefficients except the coefficient of dx in / 2 , which we call p 2 - Consider 
the set of isolated solutions over this line yields a variety of dimension 1 which we intersect with TZi 
defined by 

Ki = {p2-ri){x-r2) 

for random ri,r 2 G C. Using Bertini, we obtain 132,091 points, which was verified to be complete by 
the trace test, with exactly 79,135 satisfying p 2 — ri = 0. Therefore, our computation shows that the 
BKK (polyhedral) bound for this system is 79,135. Since this contradicts the bound of 83,977 reported 
in 1391, we provide a proof based on exact computations in polymake |10l that 79,135 is indeed correct. 

Theorem 7.1 (BKK bound for Alt's problem). The BKK (polyhedral) bound for fi,... ,/i 2 as above for Alt's 
problem is equal to 79,135. 

Proof. Following 139 1, we can use fi,..., fi to remove n, n, m, and m resulting in an unmixed system, 
i.e., all polynomials have the same monomial support, consisting of 8 polynomials of degree 7 in 8 
variables with the same BKK (polyhedral) bound as the original system. Since 7 y, and 7 ? are quartic 
polynomials, one naively would have expected the polynomials to have degree 8 , which is not the case 
due to exact cancellation in the coefficients. Hence, to properly recover the monomial support, we used 
the following computation in Maple to find that the support of these 8 polynomials is 239 monomials. 

#input 

n := a*xHat: nHat := aHat*x: m := b*yHat: mHat := bHat*y: 

px := nHat - x*dHat: py := mHat - y*dHat: 

qx := n - xHat*d: qy := m - yHat*d: 

rx := d*(aHat - xHat) + dHat*(a - x) - d*dHat: 

ry := d*(bHat - yHat) + dHat*(b - y) - d*dHat: 

g := qx*ry - qy*rx: 
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gHat := rx*py - ry*px: 

gZero := px*qy - py*qx: 

f := g*gHat + g*gZero + gHat+gZero: 

nops([coeffs(expand(f),[a,b,x,y,aHat,bHat,xHat,yHat])]); 
239 #output 


We then used the software polymake lUOl to compute the vertices of the polytope which is the convex 
hull of these 239 monomials resulting in 150 vertices. We note that 1391 reported 259 monomials and 158 
vertices with the lower values in our computation possibly resulting from using symbolic computations 
in Maple to remove monomials whose coefficients are identically zero. To complete the proof, we used 
polymake to compute the volume of this polytope, which was 2261/1152. Hence, the BKK (polyhedral) 
bound is equal to 

8! -2261/1152 = 79,135. 


□ 
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